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Summary: 

High power multiple quantum well AIGaAs diode laser master oscillator - power 
amplifier (MOPA) systems were examined both experimentally and theoretically. 
For two pass operation, it was found that powers in excess of 0.3 W per 100 pm of 
facet length were achievable while maintaining diffraction-limited beam quality. 
Internal electrical-to-optical conversion efficiencies as high as 25% were observed at 
an internal amplifier gain of 9 dB. Theoretical modeling of multiple quantum well 
amplifiers was done using appropriate rate equations and a heuristic model of the 
carrier density dependent gain. The model gave a qualitative agreement with the 
experimental results. In addition, the model allowed exploration of a wider design 
space for the amplifiers. The model predicted that internal electrical-to-optical 
conversion efficiencies in excess of 50% should be achievable with careful system 
design. The model predicted that no global optimum design exists, but gain, 
efficiency, and optical confinement (coupling efficiency) can be mutually adjusted 
to meet a specific system requirment. A three quantum well, low optical 
confinement amplifier was fabricated using molecular beam epitaxial growth. 
Coherent beam combining of two high power amplifiers injected from a common 
master oscillator was also examined. Coherent beam combining with an efficiency 
of 93% resulted in a single beam having diffraction-limited characteristics. This 
beam combining efficiency is a world record result for such a system. Interferometric 
observations of the output of the amplifier indicated that spatial mode matching 
was a signifcant factor in the less than perfect beam combining. Finally, the system 
issues of arrays of amplifiers in a coherent beam combining system were 
investigated. Based upon experimentally observed parameters coherent beam 
combining could result in a megawatt-scale coherent beam with a 10% electrical-to- 
optical conversion efficiency. 
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Ridge Guide Semiconductor Diode Laser 



p-GaAs cap ■ GaAs active layer and waveguide 

■I far-field emission 

■ p-AIGaAs cladding and ridge 
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A broad-area AlGaAs amplifier operating cw has delivered 425 mW of total power with 342 mW in a single lobe 
diverging at 1.02X the diffraction limit (FWHM 0.483°, 87,4-pm actual aperture) for a master oscillator input power 
of 70 mW. The spatial coherence of the amplifier output is 0.97, and the mutual spatial coherence between the 
oscillator and amplifier is 0.96. 


Recent reports have demonstrated considerable sin- 
gle-lobe power for master oscillator-power amplifier 
(MOPA) systems including a 160-Mm-wide injection- 
locked amplifier having M).45 W in a single lobe with a 
divergence 1.1X the diffraction limit, 1 a 400-^m-wide 
amplifier developing 0.51 W in a single lobe with a 
divergence 1.5X the diffraction limit with the use of 
two-stage injection locking, 2 and a 500-pm-wide trav- 
eling-wave amplifier having 1 W in a lobe with a diver- 
gence ~3.6X the diffraction limit. 3 Though the devi- 
ation from the diffraction limit is a useful gauge for the 
wave-front quality of the output beam, the spatial 
coherence 4 is more fundamentally related to the effi- 
ciency of wave-front transformations, the ability of 
the beam to carry information, and the ability to com- 
bine several beams coherently. It has recently been 
pointed out that the divergence is not a sensitive mea- 
sure of the coherence. 5 There are numerous reports 
on the mutual coherence of the elements of diode-laser 
arrays, 5-5 a report for injection-locked arrays, 9 and a 
report for a traveling- wave amplifier array. 10 The 
highest values of the coherence, as represented by the 
fringe visibility V have been >0.86 across a 20-element 
laser array, 0.6 for the elements of an injection -locked 
array, and 0.9 for the elements of a traveling- wave 
amplifier array. A visibility V * 0.87 has been report- 
ed for the mutual coherence between a diode-laser 
master oscillator and a single-stripe injection -locked 
amplifier. 11 

The research reported here demonstrates a 96-Mm- 
wide amplifier developing 342 mW in a single lobe that 
diverges at only 1.02X the diffraction limit, a spatial 
self-coherence for the amplifier of 0.97, and mutual 
spatial coherence between the oscillator and amplifier 
of 0.96. 

The interference of two quasi-monochromatic 
beams of intensities /i and I 2 at positions within the 
respective beams of xi and *2 provides information 
about the coherence between the two beams. The 
experimentally measurable fringe visibility V at some 
observation point P in the overlap region between the 
two beams is given by 5 


s Anax ^min ^ 

Anax Anm A(*l) ^2^2^ 

x W*i. *2. t)|, (1) 

where 712 is the complex degree of coherence for 
beams 1 and 2 that depends on the relative positions 
within the beams at which they overlap at P and the 
time delay r between the beams. For r much less than 
the coherence time, 712 represents the spatial coher- 
ence of the beams. When Ii = / 2 the fringe visibility 
gives directly the magnitude of the mutual coherence 
712. Some limiting cases of the degree of coherence 
that are of interest are the mutual coherence between 
two different beams, 

0 < l7i 2 (*i, x 2 )\ < 1, (2) 

and the self-coherence, i.e., the coherence between 
portions of the same wave front, 

= *r)l * 1, (3) 

0 ^ l7n'(*i ^ 3Cj')l — 1. (4) 

The MOPA system consisted of a single-mode 
AlGaAs diode laser (Spectra Diode Laboratories 
SDL5410) delivering 70 mW to the amplifier coupling 
optics. The double-pass amplifier (Spectra Diode 
Laboratories SDL-2419) was a 10-stripe gain-guided 
array that was 96 Mm wide and 250 Mm long. This is 
referred to as a broad-area amplifier based on findings 
that the gain-guided array behaves like a broad-area 
gain medium with periodic perturbations in the gain 
due to the localized current sources provided by the 
array stripes. 12 The rear facet had a high-reflectivity 
coating, and the front facet had an antireflection coat- 
ing. The antireflection coating on the front facet 
raised the threshold current to ~620 mA from a value 
of 287 mA for the uncoated facet. The anamorphic 
coupling optics were similar to an earlier paper. 13 The 
oscillator lateral beam shape at the entrance face of 
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Fig. 1. Far field of the amplifier with 425-mW output pow- 
er (lower trace). The strong lobe has a divergence (FWHM) 
of 0.483° , and the power contained within the lobe is 342 
mW. The lobe width is 1.02X the diffraction limit for the 
measured width of the near field, 87.4 ^m. The upper trace 
is the normalized integral of the far-field intensity showing 
80% of the total intensity in the strong lobe. 


the amplifier was approximately Gaussian with a 
width (FWHM) of 120 pm. The oscillator beam was 
incident at 5° to the amplifier axis so that the ampli- 
fied beam exited through the spherical and cylindrical 
coupling optics and could be picked off by a mirror. 
Far-field measurements of the lateral beam profile 
were made with a linear charge-coupled -device array. 
Near-field measurements of the amplifer were made 
by picking off the amplified beam in front of the cylin- 
drical lens. A Faraday isolator between the oscillator 
and amplifier provided 35 dB of isolation. The oscil- 
lator and amplifier were both operated cw with the use 
of stabilized current sources and held at constant tem- 
perature with the use of thermoelectric cooling. 

Coherence mesurements were made with a shearing 
interferometer that had equal path lengths (i.e., path- 
length difference much less than master oscillator co- 
herence length) in order to ensure that the fringe- 
visibility measurements were a function of only the 
spatial coherence. The initial setup of the interfer- 
ometer gave stable straight-line fringes perpendicular 
to the p-n junction planes. A phase-shifting method 
proved most accurate for the determination of the 
visibility. A photomultiplier was placed behind a 25- 
^m pinhole, which was then placed in the overlap 
region of the beams from the two arms of the interfer- 
ometer. The individual intensities of the two beams 
transmitted trough the pinhole were adjusted to be 
equal. A mirror in one arm of the interferometer was 


mounted on a piezoelectric stack and driven by a si- 
nusoidal voltage waveform to vary the path length of 
that arm of the interferometer by several wavelengths. 
The intensity variations as several fringes moved 
across the pinhole were recorded on a digital oscillo- 
scope. The maximum and minimum intensities re- 
corded were used in the calculation of the fringe visi- 
bility. A series of measurements were made by main- 
taining one of the beams fixed in position and scan- 
ning the beam from the other arm of the interferome- 
ter across the pinhole; visibility measurements were 
made for each setting. All measurements were made 
with fringe spacings large enough that spatial averag- 
ing did not reduce the measured visibility. 

The total output power of the amplifier was 425 m W 
with a 70-mW signal from the master oscillator into 
the amplifier coupling optics and a 671-mA current to 
the amplifier. The far-field profile for these condi- 
tions is shown in Fig. 1. The predominant single lobe 
with far-field divergence of 0.483° (FWHM) contains 
80% of the total power, or 342 mW, as is shown by the 
integral of the far field in Fig. 1. The measured near- 
field profile was 87.4 fi m in width. This is narrower 
than the 96 -jim actual current-pumped width of the 
amplifier because of vignetting by the lossy regions on 
either side of the amplifier. A lobe width of 0.483° i 9 
1.02X the Fraunhofer diffraction limit for a uniformly 
illuminated 87.4-jim slit. The other features in the far 
field are also readily understood. The small lobe at 
2.5° is the result of diffraction by the grating formed 
by the gain modulation created by the 10 stripes in the 
amplifier. 12 ’ 14 The small sidelobes at 7.7° and 8.3° are 
features that were observable in the near field of the 
master oscillator and are thus not the result of beam 
degradation due to the amplifier. The amplified 
spontaneous emission power along the amplifier axis 
was <30 mW. 

Four sets of fringe- visibility measurements were 
made to characterize the spatial coherence of our 
MOPA system. The measurements were made with 
100-mW power out of the master oscillator, 65-mW 
power incident to the amplifier, and 400-mW total 
power out of the amplifier. The fringe visibility of the 
oscillator (designated by the subscript 1) and amplifi- 
er (designated by the subscript 2) were measured 
across their beam profiles x? with respect to the beam 
center * 0 of either a split-off portion of itself 



Fig. 2. Intensity fringes obtained from the phase -shifting 
shearing interferometer. This measurement is of l 7 l2 (xi = 
0, x 2 30 0)1 and gives a value of 0.97. 




[une 15, 1991 / Vol. 16, No. 12 / OPTICS LETTERS 


915 



Fig. 3. Measured coherence as a function of the spatial 
position in one beam for the four coherences described in the 
text. Representative error bars are shown on one point. 
The data points are at the beam center and at the 75%, 50%, 
30%, and 10% intensity points. 

[Ith'(*i - o, l 722'(*2 = o, * 2 ')l] or the other beam 
[|y 2 i(x 2 - 0, xi)[, b-i 2 (*i * 0, x 2 )l]. An example of the 
experimentally acquired fringes for the measurement 
ofl 7 i 2 (*i * 0, x 2 = 0)1 is shown in Fig. 2. The fringe 
visibilities for the four sets of coherences are shown in 
Fig. 3, measured at the beam center and at the 75%, 
50%, 30%, and 10% intensity points. The position 
parameter x,- is normalized to the beam HWHM. The 
observed visibility ranges from 0.98 to 0.95. The val- 
ues of lyirfxi ■ 0, xy * 0)1 * 0.98 and l 7 22 '(x 2 * 0, x 2 ' m 
0)1 m 0.97 both differ from the ideal expectation of 1.0. 
Two experimental factors, (1) the inability to identify 
and overlap the corresponding points on the split wave 
front perfectly and (2) light scattering from dust in the 
two arms of the interferometer, are thought to be the 
causes of the deviation of the measured values from 
the ideal. In all cases there is a general decrease in the 
observed fringe visibility at the edges of the beam. A 
contribution of a few percent or less of a higher-order 
spatial mode in the master oscillator would be suffi- 
cient to explain this observation fully. 15 Other fac- 
tors, such as increased contribution of spontaneous 
emission to the intensity, might also contribute to the 
slight loss of coherence away from the center of the 
beam. A useful measure of the fraction of the light in 
one be am that is coherent with a point in the second 
beam is the intensity- weighted coherence T;/: 

Tjj = ^ = *;>l I (5) 

where Ij(xj) is the intensity of the beam at a position in 
the profile xj and Ax, is the measurement interval. 
For the four cases examined here, Tn = 0.97, r 2 2 “ 
0.97, r 12 « 0.96, and T 2 i - 0.96. The loss of 0.01 in the 
mutual coherence relative to the self -coherence seems 
surprisingly small in light of earlier results. 5 ^ 11 Possi- 
ble explanations for any loss of coherence include the 
experimental factor of light scattering from the addi- 


tional optics coupling into and out of the amplifier arm 
of the interferometer and the more fundamental co- 
herence degradations due to amplified spontaneous 
emission and the noise dynamics of the amplifier. In 
any case the small loss in coherence induced by this 
high-power amplifier indicates that multiple amplifier 
chains 2 and amplifier arrays 9 16 17 are practical tech- 
niques for generating highly coherent beams of high 
power. 

The results described here demonstrate that high- 
power diffraction-limited single-lobe output with high 
coherence can be obtained from a single-mode master 
oscillator, broad-area power amplifier system. The 
spatial-coherence values reported here are, to our 
knowledge, the highest reported for injection-locked 
or traveling-wave MOPA’s. The high power and high 
spatial coherence are particularly important for the 
realization of coherent beam-combining systems for 
power scaling. 1617 

Note added in proof: A paper published after sub- 
mission of this Letter has reported values of r 12 = 0.96 
for a pair of injection-locked amplifiers. 18 
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Coherent summation of the output of two saturated traveling-wave amplifiers, each injected from a common 
semiconductor laser and all operated cw, has resulted in a single diffraction -limited beam containing 93% of 
the power originally contained in the two individual amplified beams. The 372 mW of power in a single 
diffraction- limited lobe is a 1.5-fold enhancement over the power available from a single amplifier. 


Semiconductor laser arrays 1 and broad-area ampli- 
fiers 23 can now deliver powers in the range 1-3 W 
cw while maintaining high spatial coherence. Ad- 
vances in coherent summation are complementary 
to efforts at increasing the power of semiconductor 
lasers and amplifiers in that coherent summation 
offers a means, independent of specific device 
optimizations, to combine the power from several 
mutually coherent high-power sources into a single 
beam. Three families of coherent summation meth- 
ods have been explored: collinear interferometric 
summation, 4-8 aperture filling, 9-15 and two-wave 
mixing. 16,17 The mutual coherence among the beams 
to be summed can be created by placing all the gain 
elements inside a laser resonator 4,5,1011 or through 
extracavity amplification of mutually coherent wave 
fronts. 6-8,12-15 The three summation methods share 
the requirement of mutual coherence and wave-front 
matching of the beams to be combined. 

Among the collinear interferometric summation 
methods, binary phase gratings have been used to 
sum two injection-locked Nd:YAG lasers with an effi- 
ciency of 0.75. 6 Two-wave mixing of two injection- 
locked AlGaAs lasers in BaTi0 3 has resulted in 
summation efficiencies as high as 0.80, yielding 
98 mW of power in a diffraction-limited single 
lobe. 16 The interferometric power amplifier (IPA), 
consisting of a Mach-Zehnder interferometer with a 
semiconductor amplifier in each interferometer arm, 
has demonstrated a combining efficiency of 0.83, for 
a total output power of 8.3 mW (Ref. 7) and 110 mW 
with a combining efficiency of 0.85. 8 

In this Letter we report on experiments demon- 
strating collinear coherent summation with a Mach- 
Zehnder IPA. The IPA sums the output of two ampli- 
fiers into a diffraction-limited lobe with an efficiency 
of 0.93 and an output power in the single lobe of 
372 mW. This is a 1.5-fold power enhancement over 
the power obtainable from a single amplifier. To our 
knowledge, this is the highest-efficiency interferomet- 
ric summation with semiconductor gain elements and 
the highest single-lobe power obtained by interferom- 
etric summation reported to date. 

0146-9592/93/080619-03$5.00/0 


Collinear interferometric summation 7 can be ac- 
complished when two mutually coherent beams with 
intensities I x and / 2 intersect on the surface of a beam 
splitter (Fig. 1). The beam-splitter surface bisects 
the intersection angle between the two beams, result- 
ing in two beams I A and Ib exiting the beam splitter 
that are the coherent sum of the two input beams, 

Ia = I\R + I 2 T + 2yi2>/f 1RI2T co^A <f>) , (1) 

Ib = I\T 9- I 2 R 9- 2yi2>/7i7V^R cos( A 9- 77 ), (2) 

where R and T are the reflectivity and transmission 
of the beam splitter, yi 2 is the degree of mutual 
coherence between the two beams 1 and 2, and A<£ 
is the phase difference between the two beams. If 
A <j> - 0, yi 2 = 1, R — T f and /1 = / 2 , then 

I a - h + h , ( 3 ) 

I B = 0 . (4) 



paths is an indication of the relative intensities in each 
path. The beam splitters are designated BS1 and BS2. 
The coupling optics are not shown. 
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The mutually coherent beams 1 and 2 are combined 
into a single beam A containing all the power origi- 
nally in beams 1 and 2. If y 12 < 1 , the spatial-mode 
profiles of beams 1 and 2 are not perfectly matched, 
or the intensities and beam splitting are not properly 
matched, 7 the summation will not be complete. A 
measure of the efficiency fl of summation is 



To understand how the IPA can lead to power 
enhancements over a single amplifier, it is necessary 
to examine the power output of amplifiers with sat- 
urated gain. 9 Compare the power output of a single 
amplifier P amp with the power output of a two-arm 
IPA P IPA for identical injection power from the master 
oscillator P in and all individual amplifiers operating 
at the same small-signal gain G 0 : 

^amp = Psat ^ l n (G 0 /G), (6) 

where P 8at is the saturation power, the actual gain 
G is 

G - Pamp/An , (7) 

and P in is the injected power. For the IPA with two 
identical amplifiers and equal beam splitting, 

PiPA = 2flP amp [/TP m /2)L (8) 

where P amp [f(PJ2)] describes the use of Eq. (6), 
where PJ2 replaces P m in Eq. (7). From Eqs. (6)- 
(8), it is apparent that, when the single amplifier 
is unsaturated, it will deliver more output power 
than the IPA when the coherent summation process 
has an efficiency fl < 1. However, as the injected 
power increases, both the single amplifier and the IPA 
saturate. Because the injection power is distributed 
among the amplifiers of the IPA, the IPA will undergo 
less gain saturation. Thus the gain versus injected 
power curves for the IPA and the single amplifier 
cross at some value of the gain. This crossover 
point can be derived for the ideal case of identical 
amplifiers and equal beam splitting if Eqs. (6) and (8) 
are set equal and rearranged to provide the following 
relationship: 



For greater gain saturation than that given by Eq. (9), 
the IPA will have higher-power delivering capabilities 
than a single amplifier. 

The experimental schematic is shown in Fig. 1. 
The master oscillator was a single-mode AlGaAs 
laser, and all the amplifiers were 10 -stripe gain- 
guided arrays with an aperture width of 100 /xm 
and a length of 250 ju m. The amplifier had nominal 
front-facet reflectivity = 0.03 and rear-facet 
reflectivity R 2 = 0.97. The apparent threshold 
current was >600 mA. The master oscillator and 
all gain elements were operated cw. The amplifiers 
were operated in a double pass, 5° off axis. The 


preamplifier output was split into two beams at 
BS1. A mirror mounted on a piezoelectric stack 
in one arm of the interferometer provided relative 
phase adjustment at BS2. The outputs from the two 
amplifiers were coherently summed at the final beam 
splitter BS2 when the relative phase delay between 
the two arms was set to zero. The summed output 
at A was monitored by a line-scan CCD camera. 
The amplifiers used in these experiments, although 
nominally the same as the amplifier described in 
Ref. 18, were only capable of lower individual output 
powers for coherent operation due to individual 
device variations. 

The far field of the coherently summed output 
of amplifiers 1 and 2 is shown in Fig. 2, curve a. 
The in-phase far field had a maximum output power 
of 372 mW in a diffraction-limited lobe. The far 
field for the out-of-phase condition shows nearly com- 
plete cancellation of the intensity (Fig. 2, curve d) 
demonstrating the high quality of the mode match- 
ing and high mutual coherence of the two amplified 
beams. The far fields of the individual amplifiers 
had diffraction-limited central lobes containing 92% 
of the total power for amplifier 1 and 80% of the 
total power for amplifier 2. The far fields from the 
amplifiers in each of the arms of the interferometer 
are shown in curves b and c of Fig. 2. Before BS2, 
the amplifier output powers in the central lobe were 
206 and 195 mW for amplifier currents of 600 mA 
each. The corresponding interferometer combining 
efficiency was H = 0.93. The divergence of the in- 
phase summed lobe is equivalent to that of the in- 
dividual amplifiers. The individual amplifiers were 
operated at overall gains of 7.0 and 6.0 dB, whereas 
the small signal gains were 12.3 and 14.3 dB. Based 
on the estimated coupling efficiency, the internal 
gains are 10 dB higher than the overall gains. 

The IPA output power for the in-phase condition 
(coherent summation) is shown in Fig. 3 as a func- 
tion of input power. This can be compared with the 



Relative angle, deg 


Fig. 2. Far-field profiles for the IPA. The vertical scale 
is identical for all profiles, and the baselines are offset for 
clarity. Curve a, coherently summed far field from two 
amplifiers with 372-mW single-lobe power; curve b, far 
field of amplifier 1 with 206-mW single-lobe power; curve 
c, far field of amplifier 2 with 195-mW single-lobe power; 
curve d, coherently summed out-of-phase profile. 
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Fig. 3. Total output power of a single amplifier (O) op- 
erated at 550 mA and the output power of the IPA (□) as 
a function of the total input power. The solid curve is a 
theoretical fit for the individual amplifier that includes 
gain saturation, and the dashed curve is a theoretical 
fit for the IPA power that includes the efficiency of the 
coherent summation. 



0.01 01 1 10 100 1000 
Input Power (mW) 

Fig. 4. Gain saturation of the single amplifier (A) and 
the IPA (□). The curves are theoretical fits as in Fig. 3. 

output power of a single amplifier (before spatial 
filtering of the diffractive sidelobes) (Fig. 3). At high 
output powers, the IPA gives up to 1.3 times the 
total power and 1.5 times the single-lobe power of 
an individual amplifier for a given total injected 
power owing to reduced gain saturation. The data 
shown in Fig. 3 are replotted as gain versus the 
logarithm of the input power in Fig. 4 to demonstrate 
more clearly the crossover predicted by Eqs. (6)— (8). 
The calculated curves shown in Figs. 3 and 4 use 
the extensions of Eqs. (6)- (8) that include amplifier 
internal losses 19 and fl = 0.93. The vertical distance 
between the two saturation curves in Fig. 4 allows 
determination of the power-output enhancement as 


a function of the gam saturation. Deeper saturation 
leads to a larger power enhancement for the IPA, i.e., 
a larger difference in the gam. Because of the non- 
linear characteristics of the gain saturation, doubling 
the number of amplifiers does not necessarily lead 
to a doubling of the output power for a given input 
power. 

Coherent summation of two AlGaAs amplifiers in 
an interferometric power-amplifier configuration has 
been done with 93% summation efficiency and re- 
sulted in 372 mW of power in a single diffraction- 
limited lobe. Coherent summation of an array of 
gain elements provides a means to increase the avail- 
able optical power by at least a factor of 1.5 over a 
single equivalent gain element. 
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Symbol definitions 


Symbol 

Definition 

Vajjj£ 

Units 

gs(z) 

modal gain as a function of position 


cm'* 

N d °(z) 

transparency electron density at reference temperature 


cm*^ 

N d °(T) 

temperature dependent transparency electron density 


cm'^ 

N w 

number of quantum wells 



r 

optical confinement per well 



Go 

gain constant 


cm’* 

X s P 

spontaneous emission lifetime 


s 

d 

thickness of junction 


cm 

e 

electron charge 

1.6 10* 19 

Coulomb 

J 

current density 


A/cm" 

c 

speed of light in vacuum 

3.0 10 10 

cm/s 

n g 

group index 

3.3 


N d (z) 

electron density as a function of position 


cm'^ 

S(z) 

photon density as a function of position 


.3 

cm 

AT 

change in temperature relative to the reference value 


°K 

To 

characteristic temperature 

150 

°K 

R th 

thermal resistance 


°K/VV 

R s 

serial resistance 


Q 

Vf 

forward voltage 

1.5 

V 

L 

length of the amplifier 


cm 

D 

width of the amplifier 


cm 



Symbol 


Definition 


Mas. 


Units 


Pd 

electrical power fed to the amplifier 

W 

Pin 

optical power injected to the amplifier 

W 

P «M2 

optical power transmitted from facet 2 

w 

p outl 

optical power transmitted from facet 1 

w 

p av 

fraction of electrical power transforming to heat 

w 

p nHX 

facet damage limit for power density 

W/cm 2 

G s 

single pass gain 


Rl,R2 

front and back facet reflectivities 


P 

optical mode propagation constant 

cm'* 

h 

planck’s constant 

cm’* 

a 

distributed losses 

cm 

X 

lasing wavelength 


Jo 

transparency current density 

A/cm^ 

b 

gain - current coefficient 

cm/A 

« (±) 

injected / reflected amplitude 


b {±) (z) 

forward/backward travelling amplitudes 


c {±) 

transmitted and amplified optical amplitude 




Laser Amplifiers 


Definition 


Optical amplifiers can be divided into two different categories: Coherent optical amplifiers and in- 
coherent optical amplifiers. Coherent optical amplifiers are devices that increase the amplitude of 
an optical field while maintaining its phase. If the input optical field is monochromatic, the output 
will also be monochromatic, with the same frequency. In contrast, an amplifier that increases the 
intensity of an optical wave without preserving the phase is called an incoherent optical amplifier. 
This thesis is concerned with coherent optical amplifiers, namely quantum well stuctures. Such 
am- 
plifiers are used in many applications; examples include the amplification of weak optical pulses 
for use in long distance fiber communications, the amplification of highly intense optical pulses for 
laser-fusion applications or the amplification of continous wave optical field from a master oscil- 
lator in a MOPA system. 

Coherent amplification is achieved through amplification by the stimulated emission of radiation. 
That is a photon in a given mode induce an atom at a higher energy level to undergo a transition to 
a lower energy level, emitting a second photon into the same mode (a photon with the same direc- 
tion,frequency and polarization). The two photons, in turn, can stimulate the emission of two ad- 
ditional photons with the same properties and so on. The requirement for stimulated emission being 
the incident photon to have a nearly equal energy to the corresponding atomic transition, the pro- 
cess is restricted to a band of frequencies determined by the atomic linewidth. Thus, the primary 
frequency selection in optical amplifiers is due to the energy level differences of the laser material 



used. Optical cavities may be used to assure auxiliary frequency tuning in contrary o electronic 
amplifiers where the frequency tuning is achieved through an internal resonant circuit making use 
of inductors and capacitors. 

In thermal equilibrium the number of atoms in the lower energy level is higher than the number of 
atoms in a higher level, resulting into the domination of absorption by the low energy atoms over 
the stimulated emmision through the minority high energy atoms. To achieve amplification, the 
inequality mentioned above has to be reversed with the help of an external source of power result- 
ingfinto a nonequilibrium situation. An elementary diagram describing the operation of a laser am- 
plifier is illustrated below (see figure 1.1). 


PUMP 


ATOMS 



Figure 1.1 : Simplified illustration of amplification along a traveling wave amplifier 


Caracterization criterias for optical amplifiers 

Real coherent amplifiers are caracterized by a gain and phase shift that are frequency dependent. 
Furthermore, for a sufficiently large input photon density the amplifier exhibit saturation which re- 
sults into a non-proportional relationship between input and output powers. Also saturation may 



result int<i generation of harmonic components. Finally a noise is always present in the output. 
An amplifier may therefore be characterized by the fallowing: 


- Gain 

- Bandwidth 

- Phase shift 

- Power source 

- Nonlinearity and gain saturation 

- Noise 


Examples of laser amplifiers 

Laser amplification can happen in a variety of materials. Usually each system has interacting ener- 
gy levels influencing the electron populations of the transition of interest. The operations of most 
of these amplifiers can be modelled as a three or a four level system. Examples are the three-level 
ruby laser amplifier, the four-level neodymium-doped yttrium-aliininium garnet laser amplifier, 
and the three-level erbium-doped silica fiber laser amplifier. Besides the solid state lasers men- 
tioned above gas and liquid lasers can also be used to achieve amplification. 

Pumping may be implemented using many methods, including the use of electrical, optical and 
chemical means. Gas lasers usually requires a direct current (dc) or a radio-frequency (RF) dis- 
charge current for pumping. The solid-state silica fiber laser requires optical pumping which may 
be achieved through a semiconductor' in lection diodes Finally semiconductor laser amplifiers use 
direct current (dc) or ftlfarnativetcurreni (acV’for pumping in a cw or pulsed mode of operation res- 



pectively. 


Optical Amplification bv semiconductor lasers 

After the invention of semiconductor laser [1] in 1962, the coherent light properties of the device 
became the focus of extensive research [2]. Due to the inherent nature of optical amplification, 
semiconductor lasers can be utilized as optical amplifiers besides oscillators. It acts as a “ linear” 
amplifier when the supplied current is below the oscillation threshold and as a nonlinear amplifier 
through injection locking when it is operated above threshold. 

In the linear mode, it is possible to classify the amplifier into two types: the Fabry-Perot (FP) and 
travelling wave (TW) type, differing from each other at the reflectivities at both end mirrors. The 
travelling wave amplifier has very low facet reflectivities resultingjmq the incident light being am- 
plified in a single pass without suffering any reflection This kind of amplified 'causes a smooth 
increase in spontaneous emission with increasing electrical power, making difficult to differentiate 
amplification below the threshold from amplification above the threshold babry-Perot type ampli- 
fiers/howcver,\are caracterized by considerable facet reflectivity resulting into a regenerative reso- 
nant amplification between both end mirrors. 

The principle underlying the operation of a semiconductor laser amplifier is the same as that for 
other laser amplifiers: population inversion. It is achieved through artificial} pumping of electrons 
to higher energy levels by electrical current injected in a p-n junction diode. Application of a for- 
ward bias voltage results into') the injection of carriers pairs into the active region , where they re- 
combine through stimulated emission. 

The theory behind semiconductor laser amplifier is somewhat more complicated than the amplifier 



structure mentioned earlier because the transitions occur between bands of closely spaced energy 
levels rather than well-separeted discrete levels. However, the semiconductor laser amplifier may 
be viewed as four level laser system (see the appendix A), in which the upper two levels lie in the 
conduction band and the lower two levels lay in the valence band. 

These amplifiers have both advantages and disadvantages. Although they are very small in size so 
that they can easily be incorporated into optoelectronic integrated circuits and their bandwidths can 
be as large as 10 TH? (greater than optical fibers), the coupling loss is very high (3dB to 5dB) and 
they are temperature and polarization sensitive. 

Motivation for Heterostructures 


The complex dependence of the gain coefficient on the injected-carrier concentration makes the 
analysis of the semiconductor amplifier somewhat difficult. Because of this, it is customary to 
adopt an empirical approach in which the peak gain coefficient g s is assumed to be linearly related 
to the injected carrier concentration N e j for values of Ngj close to the operating point (similar to the 
small signal treatment for electronic amplifiers). The dependence of the peak gain coefficient g s on 
Nfi may then be modeled by the linear equation 


g = G 

6 S o 


(N 


el 


- 1 


V*;, 


(i) 


The parameters a and 



are chosen to satisfy the requirement that when the injected carrier con- 


centration is zero, the peak gain coefficient is equal to minus aloha which is the absorption coeffi- 



cient of the semiconductor and when the injected carrier density is equal tc\ transparency current 
density, the medium becomes transparent, that is gs is equal to zero. It is clear from the equation 
that, depending upon whether we are operating below or above transparency, the amplifier acts as 
a strong absorber and strong amplifier respectively. 
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Figure 1.2 : The energy-band diagram of a forward-biased heavily doped p-n junction 

The energy-band diagram of a forward-biased heavily doped p-n junction is illustrated below (see 
figure 1 .2). The conduction and valence bands quasi-Fermi levels Ef c and Ejy lie within the con- 
duction and valence bands respectively and a state of quasi-equilibrium exists in the active region. 
The quasi-fermi 16o>els are sufficiently well separated so that a population inversion is achieved 
and net gain may be obtained over the bandwith E g <hv < Ej c - E fv within the active region. The 
thickness d of the ative region in a semiconductor homostructure laser is a very important param- 
eter that is determined principally by the diffusion lengths of the minority carriers at both sides of 



the junction. Due to the proportionality between electron densities and injected current density at 
low input optical powers, equation (1) can be rewritten as 




(2) 


where J and Jo can be expressed as 





( 3 ) 


It is seen from the above equations that transparency current density Jo is directly proportional to 
the junction thickness d, so that a narrower active region would result in a lower current density to 
achieve transparency. The amplifier however would still experience the same gain at a much lower 
current density. The new device would consequently require a much easily achievable electrical 
power and would be much more efficient. However, reducing the size of the junction causes some 
problems, namely the diffusion of electrons and holes out of this smaller region, their diffusion 
lengths being greater. A solution to the problem of confining electrons and photons into an active 
region that is smaller than their diffusion lengths is to use a heterojunction device which would also 
confine the optical beam to an active region which is smaller than its wavelength. 


Heterostructures 


The concept of the heterostructure is to form heterojunction potential bariers on both sides of the 
p-n junction to provide a potential well that limits the distance over which the minority carriers may 



diffuse. This confinement results in active regions of thickness as small as 0.1 |im ( even thinner 
confinement regions, = 0.01 were achieved with quantum wells). 

It is also possible to achieve optical confinement of the electromagnetic beam simultaneously, by 
selecting the material of the active region in such a way that its refractive index is slightly greater 
than that of the surrounding layers so that the structure acts as an optical waveguide like core gui- 
ded fiber optic cables. 

A typical heterostructure device would have three layers of different lattice-matched materials: 
Layer 1 : p-type, energy gap Egp refractive index nj 
Layer 2 : p-type, energy gap Eg 2 , refractive index n2 
Layer 3 : n-type, energy gap Eg 3 , refractive index n 3 

The energy gap of the active region is chosen to be smaller than the sandwiching layers to insure 
carrier confinement and its refractive index is chosen to be slightly higher to achieve light confine- 
ment. 

The advantages of using a heterostructure device instead of a homostructure are: 

- Increased amplifier gain for a given electrical power due to the ease of reaching transparency. 

- Increased amplifier gain due to the increase of the proportion of light being amplified resulting, 
from the enhanced optical confinementt 

- Reduced absorbtion from the cladding layers, their gap energy being larger than photon enerThe 
energy band diagram and refractive index as a function of position for a double heterostructure 
semiconductor laser is illustrated in figure 1.3. 

When the thickness of the active layer is made sufficiently small (i.e. smaller than the de Broglie 
wavelength of a thermalized electron, = 50 nm in GaAS), the understanding of the physics of the 
device necessitates quantum mechanics due to the shrinking in size. Since the energy gap of the 



surrounding layers are larger than the active region, the latter start acting as a quantum well result- 
ing unto'B different energy- momentum relationship than the usual heterostructure. 
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Figure 1.3 : The energy band diagram and refractive index as a function of position for a do 
heterostructure semiconductor laser. 




Quantum well laser amplifiers 

Except for the thichness of its active layer, quantum well lasers are identical to the conventional 
heterostructure lasers. The diagram of a graded index separate confinement GaAlAs heterojunction 
quantum well laser is illustrated in figure 1.4. 
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Figure 1.4 : Anatomy of a graded index separate confinement GaAlAs heterojunction 
quantum well laser 


The function of the graded region is to enhance the advantages of using a heterostructure with res- 
pect to a homostructure laser, that is, it enhances dielectric waveguide properties of the laser and 



confines and guides electrons and holes to the well. 

It has been mentioned earlier that quantum well lasers have a different energy-momentum relation 
ship than the bulk material. The energy levels Eg of an electron of mass m c or a hole of mass m v 
confined to a one dimensional infinite rectangular well of width d can be found by solving the 
Schro dinger equation as 




(4) 


where q= 1,2,.. and m designates an electron’s or a hole’s mass. Remark that increasing the width 
of the well results in an increase in the separation between adjacent energy levels. 

As shown in the figure 1.5, electrons and holes are confined to the active region in x direction. 
However in the y-z plane , they extend to a much larger dimension and can be treated as if they 
were in bulk semiconductor. 

The energy-momentum relationship in the conduction band can be written as 


E 


, 2.2 , 2.2 , 2.2 

h k x h k 2 h k 3 

+ 2m + 2m + 2m 
c c c 


(5) 


where kj=qjn/dj k2=q2! R l t l2, ^ 3 = 4311^3 and q 14243=1,2, Remark that since the active layer 

thickness is much smaller than the surrounding layers thicknesses, the separation between the dis- 
crete values of k in x-z plane is much more important than the ones in x-y and y-z planes. Eventu- 
ally, * 2 and ^3 ma y be approximated as a continuum in which case equation (5) can be rewritten as 




Figure 1.5 : (a) Geometry of quantum-well structure (b) Energy level Diagrams in the 
quantum well (c) Cross section of E-k relation in the direction of k 2 and k 3 


where qj=l,2,3,..., and k = + |fc 3 | . In the above equation qj is the quantum number denoting 

the subband with minimum energy E c +Egj The valence band could be treated similarly. 

For the bulk semiconductor the energy is a parabolic function of the three dimensional momentum 
vector and the difference between each value of the momentum in each direction has comparable 
magnitudes unlike the quantum well structure for which the momentum vector Jfc/takes well sepa- 
rated values. 

For the bulk semiconductor the density of states can be written as (see appendix for derivation) 

P (£) = -3 

n 


(7) 




From the equation above, the density of electrons and holes near the conduction and valence bands 
respectively, fallow) as (see appendix for derivation) 


3/2 


( 2m c y 1/2 

P c( £ ) = — 2,3 (E-E c ) W1 

Z7t it 
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for E>E C and 


3/2 


( 2 m ,y i /2 

• — i L i (E,-E) 

2 71 h 


(9) 


for E<E C , where me and mv stands for the effective averaged masses of electrons and holes res- 
pectively. 

For the case of quantum well laser, the density of state functions are obtained from the magnitude 
of the two dimensional wavevector (k2, k3) and therefore they are linearly dependent on it. The ex- 
pression for the density of state per unit volume fallows as : 
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The electron density of states as a function of the energy p c (E) is related to the density of state 
function with the relation (see appendix for details). 
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Taking the derivative of the right hand side of equation (6) with respect to the two dimensional 
wavevector k yields 


dE _ h^k (12) 

dk m c 


Putting ( 1 2) into (11) yields 
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for & > E c + and 0 otherwise for qj = 1,2,... It can be seen from the above equation that when 
E > E c + E q j the density of state per unit volume is constant for each quantum number q j. Since 
to obtain the overall density of state function one has to add the density of state functions corres- 
ponding to each quantum level qj, the former exibits a staircase like shape as shown in figure 1.6. 
This property is also true for the density of states in the valence band. 
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Figure 1.6 : Density of states for a quantum-well structure and for a bulk semiconductor 


In a quantum well laser, photons interact with electrons and photons somewhat similarly to the 
bulk semiconductor, lence bands. The transitions conserve energy and momentum between con- 
duction and valence bands. Different from the bulk semiconductor, the transitions in the quantum 
well laser also conserves the quantum number q, that it they are allowed within the band qj=i in 
the valence band and qj=i in the conduction band, where i denotes the subband number. The ex- 
pressions derived fo-gain coefficient and transition probabilities in the bulk material also applied 
to the quantum well when the bandgap energy -Eg is replaced with the energy gap between the sub- 
bands, E gq = E g + E q + E q ' and a constant density of states is used rather than one that varies as 
the square root of the energy. The total gain coefficient is obtained by summing up all subband gain 
coefficients. 

At this point we can rewrite the joint density of states function in terms of frequency as 
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when h\)>E g +E q + E' , and 0 otherwise. The transition from equation ( 1 1) is straighforward 
using the relationship 
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when the transitions from all subbands are considered the joints density of states looks like stair- 
case distribution as it is illustrated in the preceeding figure. 



Gain coefficient 


As in the case of bulk semiconductor (see Appendix B), the gain coefficient of the quantum well 
structure is given by the expression 
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where fg(x > ) is the Fermi inversion factor identical for bulk and quantum well lasers. Thus the only 
difference between the two gain coefficients is the joint density of states function. The effect of dif- 
ferent density of state function on the gain coefficient is illustrated in figure 1.7, where it is as- 
sumed that the higher levels of the laser are not excited. The step wise density of state function of 
the quantum well sructure results into a smaller peak gain and narrower bandwidth. 

The dependence of the peak gain on current density supplied to the quantum well can be infered 
from the graph. An increase in the electron density will result into the generation of excess halls 
and electrons which will increase the gap between the quasi-fermi levels. When the value of the 
current density J is high enough so that Efc-Ep barely exceeds the gap Eg] between the q=l levels, 
the medium starts experiencing gain which will sharply increase up to its saturation value given by 
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If J is further increased, the spectral width will be broadened leaving the peak gain unchanged. In 
case the amount of pumping is at such a level that the difference EjfrEjy cross over Egp, the levels 
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Figure 1.7 : Density of states, Fermi inversion factor, and gain coefficient in quantum 
well and bulk structures. 

q=2 become^ excited and the peak gain coefficient experiences a second jumjxJThe described va- 
riation of the peak gain coefficient with current density is shown in figure 1 .8. Note that if q=l lev- 
els are excited only the relationship could be considered as pseudo- logarithmic. 







The semi-classical theory presented above displays the major differences between quantum well 
and double heterostructure devices. A more elaborate treatment can be found in papers where a 
pure quantum mechanical approach is used. But, the possible (and adequate) representation of the 
behaviour of the quantum well material gain with injected current density by a semi-logarithmic- 
function makes such an approach unnecessarily complexe and timewise expensive 



Figure 1.8 : Schematic relationship between peak gain coefficient and current density 



Amplifier Model 

The amplifier model chosen has to be able to predict closely the gain and efficiency of the amplifier 
in any operating condition necessitated by the optimization algorithm without increasing the ma- 
thematical complexity with unnecessary details. The model parameters should be easily deter- 
mined from experiment and the set of equations used shouldn’t become transandental at any oper- 
ating condition. 

Rate Equations 

Mathematically, semiconductor laser amplifiers are usually described by time dependent rate equa- 
tions dealing with temporal changes in electron and photon densities in the active region of the de- 
vice ignoring their spatial variation along the amplifier. Such method is used to describe amplifiers 
with high facet reflectivities for which the spatial distribution of electrons and photons are relative- 
ly even. However, to deal adequately with any kind of reflectivity (especially with zero reflectivity 
due to anti-reflective coating at facets), one has to use time and space dependent travelling-wave 
equations. 

Rate Equations in the a bscence of excitation 

In the abscence of input optical field the dynamics of amplification can be described by the fallo- 
wing two rate equations corresponding to the variation of spontaneously emitted photon and elec- 
tron densities along the amplifier. 
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The second term on the right side of equation (18) can be related to the loss coefficient of the guided 
mode as 
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The first term in the right hand side of equation (17) is the contribution of pumping electrons to the 
electron population. The second term accounts for the rate at which the excess electron population 
decays through spontaneous emission. The last term represents the saturation effect resulting from 
the depletion of electron population fallowing the electron hole combination. 

The first term at the right hand side of equation (18) accounts for the increase in the amount of 
spontaneously emitted photons at frequency u fallowing the spontaneous decay of the correspond- 
ing proportion of electrons in the conduction band. The distribution function D v is used to assure 
that the electron distribution at all wavelengths tends to increase or decrease at all wavelengths 
when photons are emitted or absorbed (homogeneous broadening). 

The steady state solution of the rate equations can be obtained by letting the left hand sides of equa- 
tions ( 17) and (18) equal to zero. Then the electron density can be obtained from equation ( 17) as 
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The summation in the right hand side of equation (20) is multiplied by two to account for noise 
photons with two possible polarizations. 

The solution for equation (18) necessitates that the spontaneously emitted photons be expressed as 



a superposition of photon densities moving in positive and negative z-directions as 
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and that the total derivative at the left of equation ( 1 8) be decomposed as 
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This way equation (18) can be written as two differential equations, one for each propagation direc- 
tion. At steady state, the derived rate equations fallow as 
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Proper integration of equations (23) and (24) will yield the solution in the form of integral equa- 
tions. The detailed development of the solution is not essential for this paper because in the range, 
the simulations are performed, spontaneously emitted photon densities are negligeable compared 
to the injected signal density. The interested reader can find the solutions in Marcuse’s paper. 



Previously, the spontaneous emission was described as an incoherent power wave because it dis- 
tributes itself continuously over a relatively wide band of wavelengths with random phases bet- 
ween adjacent wavelength components. However the signal wave has to be treated coherently 
since its transmission is affected by resonance effects especially when the reflectivities at the fa- 
cets distorts the constructive interference within succeedingly reflected waves. 

The figure below is the simplified drawing of an amplifier of length L and facet reflectivities Rj 
and R2- The amplitudes denoted with the superscript (+) and (-) represents waves propagating in 



(+z) and (-z) directions respectively. 

The differential equations describing the variation of the signal wave along the amplifier along (+) 
and (-) z-directions are given by the fallowing travelling-wave equations. 
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Where p is the effective wavevector described as 
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The first term in the right hand side of equations (25) and (26) accounts for the change in phase 
of the optical signal as it travels along the amplifier. The second terms denote the net amplification 
of the signal amplitude. Remark that both waves will grove in amplitude in the direction of their 
superscript if the modal gain is higher than scattering losses. 

The solutions of equations (25) and (26) are 

b (+) (z) = b i+) (0) exp (-/pz) exp jjJ [$ f (x) - a] dx 
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b ( \z) = b ( ] {L)exp [-ip (L -z) ] ex/> j^J [g s (x) -a]dxj ( 30 ) 



The depletion of electron density is now due to the amplification of signal photons rather than spon- 
taneous emission. The corresponding change is reflected into the electron density as 
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Boundary conditions for the Fabrv-Perot Laser Amplifier (FPLAf 


The appropriate boundary conditions for the FPLA must take into account the reflection and trans- 
mission at the facets. At z=0, the field reflectance is taken as rj and the tranmittance as tp At z=L, 
the equivalent quantities are denoted as r% t£ The relationship between mentioned reflectances and 
transmittances to the reflectivities Rj, R 2 are established as 
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The boundary conditions at z=0 and z=L (see figure) are given by 
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c (+) = t 2 b i+) (L) 

(39) 

Using equations (28), (29), (35) and (36), one can calculate the actual value of photon amplitudes 
at the boundaries as a function of one pass gain G s , facet reflectivities Rj and R 2 , input optical 
amplitude a (+> and the phase shift due to resonance effect as 
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where the one-pass gain G s is given by 


G s = J [g s (x) - a] rfx (42) 

o 

Now, one can calculate the photon densities inside and outside of the amplifier (see figure) by ma- 
king use of equations (28), (29), (37), and (38). A better understanding of equation (39) and (40) 
would be established if one obtains them through infinite summations of reflected or transmitted 
waves at the device boundaries. 

The previously mentioned method constructs the steady-state photon amplitudes at z=0 and z=L 
by adding the resulting contribution of each reflection as shown in the figure below where b[/ + \0), 
b n ^(L), a n '~ ) and c n ^ stand for 
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Thus the effective boundary photon densities are given by the previously mentioned summations 
of individual reflections at the corresponding facet as 
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Thus the effective boundary photon densities are given by the previously mentioned summations 
of individual reflections at the corresponding facet as 
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The term in the right hand side of equation (46) is nothing but a summation from zero to infinity 
of a geometrical series of radius (r j r 2 G s exp (— 2i{3L) ) . The summation converges eventually to 
equation (39) when the one pass-gain G s is smaller than the inverse of the product of reflectances 
that is when G s < , namely when the amplifier is operating below threshold. Thus equations 

(39) and (40) are mathematically correct only when the amplifier is not lasing. 

The advantage of using the approach of regrouping all the reflected waves at one boundary into one 
expression is to simplify the problem into one which only deals with two optical densities travel- 
ling in opposite direction (see figure) without dealing with the summation of infinite partial sum- 
mations. 


Relation between po wer and photon densities 


The electron and photon densities used in the model equations are implicit model parameters des- 
cribing the amplifier behaviour at microscopic level. However these parameters are not obser- 
vable and have to be linked with their corresponding counter parts at macroscopic level. Such a 



relationship has been established between the electron density and injected current density within 
the equation (30). A similar approach can be used to relate the injected optical power to the in- 
jected photon density and to the photon densities inside the amplifier through equations (28), (29), 
(39), (40). The mentioned equations are 
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Equations (48) and (49) are easily derived using difference methods. Namely, the differential 
energy A E passing through the aperture to the active region of the amplifier during a differential 
time At assuming % 100 coupling efficiency is as it is evident from the figure 



where hv is the energy of one photon and Vgis the group velocity. Writing the energy of a photon 
in terms of its wavelength X, the group velocity in terms of the group index n g and letting At ap- 
proach to zero, one obtains equation (48). 




Calculation of amplifier efficiency 


MM 



In the model, only one-sided efficiency is considered, that is the extracted power is represented by 
(Pout2-Pin) or (Pout 1 -Pin). In the former case, the reflected power Pout 1 is a loss factor decreasing 
one sided amplifier efficiency. A practical way to minimize the reflected power is to keep facet ref- 
lectivities low, that is using the device as a one pass amplifier. Another way of eliminating effi- 
ciency loss through reflected power is to use a double pass amplifier with a very high reflectivity 
at the second facet ( = 1 ) and a very low one at the first ( = 0 ), the extracted power being (Pout 1 - 
Pin). 

The two amplifier described above differ in the value of one pass gain, the latter showing lower 
gain due increased saturation effect.The two different definitions of efficiency are presented below 
for one pass and double pass amplifiers respectively. 
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The inconvenience with the double pass amplifier is that the threshold of oscillation is reached 
more easily due to the high reflectivity at the second facet changing the amplifier into an injection 
locked laser. Also, frequency detuning effects at high pump power which are undesirable for a 
MOPA array, are more profound compared to travelling amplifier case due to the decrease in band 
width. 

However it has been reported [ 1 ] that the use of a double-pass structure results in efficiency im- 
provement in quantum well laser amplifiers. 

Gain in Multi-Quantum Well Lasers 


The exact quantum mechanical calculation of the gain in a quantum well laser is cumbersome and 
computing wise very expensive. However, it is possible to use simpler expressions like the one gi- 
ven in equation ( 1 ) representing the gain-current relationship in a semi-logarithmic way. Another 
approximation developed by Wilcox [2] et al. reported good agreement between semi-logarithmic 
and calculated curves (rms error less than 5%). 
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where b and j Q have meanings of gain-current coefficient and current per well at transparency re- 



spectively. Corresponding plots illustrating the equations (1) and (20) are given below. 


GAIN IN MULTIPLE QUANTUM WELLS 



It can be noted that the two graphs resulting from two different approximation schemes have non- 
negligeable differences between them. The main reason of this difference is that the basic assump- 
tions on which quantum mechanical analysis is based are different in both cases as well as some 
structural parameters such as the aluminum content in and out of the quantum well region. Al- 
though such a variation won’t have a big impact on the theoretical optimization presented in this 
report, it may play a very important role in trying to validate the theoretical data with experimental 
one. Correcting the semi-logarithmic material gain equation with the help of some experiments will 
be necessary before starting numerical simulations.The numerical analysis presented in this paper 
makes use of equation ( 1 ) only. 

It is also important to note the dependence of transparency electron density on temperature in both 
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(1) and (22). The relationship is eventually modeled as 
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where AT is the relative temperature variation defined as 
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The serial and thermal resistances in (4) can be expressed in terms of a bulk component inversely 
proportional to the area of the diode and a constant term due to metallic contacts as 
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Equation (4) basically states that the energy not used in the amplification of the optical signal is 
converted in phonons. 

In an actual amplifier many parameters are effected by changes in temperature. However, it is pos- 
sible to lump the effect in one parameter, namely the transparency current density with the use of 
a characteristic temperature Tq. The validity of the approach for accurate modeling will hopefully 
be evaluated with the results of this thesis. However, since in semiconductors phonon recombina- 
tion results in a small number of minority carriers (unlike superconductors) the approach is com- 
monly used in the literature. 



Numerical implementation of the amplifier model 


The system of non-linear equations describing the laser amplifier model can not be solved analy- 
tically due to its mathematical complexity. The only solution method resides in numerical analysis. 
From the physics of the device inherent within the system of equations, one can assume that below 
threshold, there is a unique solution for each of the unknown variables. Above threshold, the sys- 
tem of equations becomes transandental, however that mode of operation (lasing) is out of the 
scope of this thesis. 

The aim of the simulation is to obtain signal photon densities along the amplifier lenght in order to 
determine the output power, consequently the net gain and efficiency of the amplifier. However, 
this also necessitates the determination of the variation of electron density and modal gain along 
the device together with the overall efficiency to take the effect of temperature shift into account. 
In the discrete setup, the amplifier is divided into a grid of 20 equidistant points (this sampling 
seemed to be adequate even for amplifiers as long as 0.15 milimeter) and the integrations of the 
modal gain is achieved through the trapazeoidal rule at the corresponding points of the grid. 

In the most general case, the main variables involved in the analysis are electron density (30), pho- 
ton densities corresponding to (absolute values) (28) and (29), absolute value of boundary values 
of signal amplitudes (39) (40), the efficiency of the amplifier defined as in (16) or (17), the temper- 
ature shift, and the modal gain of the amplifier ( 1 ). Among these variables, the electron and photon 
densities together with modal gain can be considered as vectors of size equal to the number of grid 
points (20). Thus, a solution vector to the problem should satisfy self-consistency at every point of 
the grid. 

The numerical analysis established is better understood if it is explained for the simplest amplifier 
structure ( travelling-wave amplifier ) first in the abscence of effects created by varying tempera- 



ture. Then, it is easier to fallow the way with which convergence problems arising from tempera- 
ture shift, high saturation, facet reflectivities and operation near threshold are handled. 


Travelling wave type amplifiers 


a)Small signal approach 


J 


l t t l i 



z 




In the case of a travelling wave amplifier, facet reflectivities are equal to zero. Consequently, the 
incident photon signal suffers no internal or external reflections and is amplified during its unique 
pass along the device The boundary photon density values within the amplifieras described in (39) 
and (40) then becomes 
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Obviously, in this structure, the total gain is equal to one-pass gain of the amplifier and there is no 
signal wave propagating in (-z) direction. The inputs to the diode are the current density J and the 
input signal power P in . The known device parameters are the length L, the width W, the number 
of quantum wells nw, the confinement factor per well T, the optical loss coefficient a, the signal 
wavelength X, the group index n g , the gain coefficient G 0 , and transparency current density J Q . 
The basic unknowns are electron and photon densities at each grid point in the amplifier. 

One can easily see in equations (30)(28) and (1) that electron and photon densities within an am- 
plifier are strongly coupled to each other. In the small signal approximation, the unknown electron 
density vector is close to its unsaturated value, making the modal gain vector close to its actual val- 
ue as well. Thus, it is possible to force the electron and photon densities to converge to a self-con- 
sistent solution by starting from the initial value of the electron density and calculating the modal 
gain and photon densities iteratively at each point of the grid by plugging one variable in the ex- 
pression of the other according to (1),(28) and (30). The method of solution is algorithmically de- 
scribed below. 

It is important to note that, self-consistency of the solution is based on the relative variation of elec- 
tron and photon densities with respect to their value in the previous iteration. If the relative varia- 
tion for each variable at every grid point is smaller than a preset infinitesimal limit (10'^), then 
convergence is achieved. It is possible to push this limit further down, however the precision of 
digital calculations and time constraints impose a lower limit. 

The described algorithm work very fast because the number of iterations needed is small due to the 
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method used and to the initial value assigned to the unknowns being close to their actual value. 
Due to small-signal approximation, the second expression in the right side of equation (30) is 
smaller than the first one, keeping the electron density positive. However, for a moderately high 
product of modal gain and photon density at given grid point, the electron density may become 
negative. This is undesirable because the the logarithm in gain expression described in ( 1 ) becomes 
undefined. To account for this problem, it is necessary to modify the convergence procedure to 
keep electron density always positive. 

b)Large signal approach 

A solution to keep electron density positive during iterations is to control its variation instead of 
letting equation (30) decides on it. The previous approach is used till the value of simulated elec- 
tron density at any point in the amplifier becomes negative. Then the electron density at that point 
is reset to its previous non-negative value. A dummy variable will also be set to the same value. At 
this point, the electron density is decreased (could be as well increased) by some percentage and 
the new value is put back in the iteration. From the calculated new gain and photon density, the 
dummy variable is calculated according to (30), the dummy variable playing the role of electron 
density. Then the value of the electron density is compared to the dummy variable. If the former is 
larger, than it is decreased by the same percentage. If it is smaller, it is increased by a lower per- 
centage and the loop is repeated again. As the iteration progresses the value of the electron density 
will approach to its actual value given by (30) without ever becoming negative. The method estab- 
lished for self-consistency is still valid with the new method. That is the electron density is accep- 
ted to converge if the percentage change in its value from one step to the other is small enough. 



Temperature effects 


It has been mentioned before that the temperature effect is lumped into one variable, namely trans- 
parency electron density. According to (55), (56) and (1), a decrease in the amplifier efficiency re- 
sults into an increase in the temperature causing the modal gain at every grid point (overall gain as 
well) to decrease reducing further the efficiency of the amplifier. 

A problem arizes while trying to incorporate this variable in the model : the efficiency of the am- 
plifier is not known until the output photon density is known. This however requires the knowledge 
of the shift in the temperature. The problem is solved iteratively assuming zero temperature shift 
in the first pass (by pass it is meant the set of model equations converging at every grid point). Then, 
the intermediate efficiency and transparency electron density of the amplifier is calculated (which 
is higher than the actual efficiency). The previously obtained electron and photon densities are 
stored in a buffer to used as initial grid values fc^the iterative steps needed to deal with the effect 
of this new variable. 

As seen in the figure, model equations are forced to converge iteratively for this new transparency 
current density starting from the initial values previously stored in the buffer. A new set of photon 
and electron density at every grid point together with a new output power and amplifier efficiency 
are calculated. Relative changes in electron and photon densities at every point and in temperature 
are calculated. Once again convergence is assumed to be achieved if relative changes are smaller 
than a predetermined value, that is the electron and photon densities in the amplifier together with 
the temperature shift don’t change significantly from one pass to the other. 

To make the model and its numerical analysis clearer, the simulation results for a travelling wave 
amplifier is plotted in figure () in terms of the variation of electron density, photon density and point 
modal gain along the amplifier. The amplifier is input an optical signal of 0. 1 5 watts and excited 
with a current density equal to 2000 amperes per centimeter square. The length and width are taken 




equal to 0.05 and 0.01 centimeters respectively. The number of quantum wells is fixed at 4, the con- 


finement factor at 1 .6 percent, the thickness of each quantum well at 75 angstrom and optical loss 














coefficient at 5 per centimeter. The expression given in (1) is used for the modal gain with gain 
coefficient equal to 1195 per centimeter and transparency current density equal to 180 amperes per 
centimeter square. 

Figure ()-i shows that the optical signal intensity grows linearly with distance within the amplifier. 
This obsevation is rather contradictory to what one might expect from equation (29) which de- 
scribes an exponential profile. The linearity can be explained with the help of the second expression 
in the right hand side of equation (31) accounting for the depletion. That is , for the amplifier sim- 
ulated, the electron density is depleted to such an extend to reduce the gain and thus the integral in 
the right hand side of (29)causing the observed linear profile. It can be shown with simple algebra 
that such linearity requires the modal point gain to vary as - as shown in figure ()-iii. The electron 
density profile fallows from the modal gain using (1). 

The analysis of the travelling wave amplifier is simplified by the abscence of optical reflections due 
to zero facet reflectivities. However, for Fabry-Perot type amplifiers, the dependence of boundary 
value signal densities on one-pass gain as described by (40) and (41) makes the analysis and nu- 
merical computations rather involved. 










Fabrv-Perot type amplifiers 

The problems encountered in the numerical analysis of Fabry-Perot type amplifiers are two-fold. 
The first difficulty is related to the presence of a signal wave travelling in -z direction and to the 
imprecision in boundary value photon densities. The second difficulty resides in the denominator 
of (40) and (41) becoming negative during iterations causing convergence problems. 

a) Non-zero facet reflectivity effects 

To clarify the way in which the effect of reflectivities is included in the model, it is assumed that 
the calculated amplifier one-pass gain resides below its theoretical threshold of oscillation value 
during computation. The way this restriction is relaxed, is described in the subsequent section. 

The methodology is somehow similar to the one adapted for temperature effects. The one-pass 
gain will initially be assumed equal to its unsaturated value. The intermediate boundary value pho- 
ton densities, given by (40) and (41), are then calculated using the intermediate one pass gain. Con- 
sequently the value of calculated photon densities within the amplifier are higher than their actual 
value. In the second pass, the one pass-gain calculated from the previous step is used in (40) and 
(41). Remark that it is now smaller than its actual value due to excessive saturation. After some 
sufficient number of passes, self consistency in electron and photon densities (thus one pass gain) 
and temperature shift is achieved. 

b) Computational problems associated with near threshold of oscillation operation 

As mentioned earlier, the system of equations describing the operation of a laser amplifier becomes 
transandental above the threshold of oscillation. During iterations, the value of one pass gain might 
be higher than its maximum theoretical value at oscillation resulting in convergence of equations 
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to a wrong set of solutions or not convergence at all. This problem is taken care similarly to the 
problem associated with negative electron density. For an amplifier for which, the denominator of 
(40) becomes negative, the variation of the denominator is no more dependent upon the one-pass 
gain. Rather, it is varied externally till it convergences to the previous denominater expression 
while remaining positive. The details of the algorithm can be seen in figure () above. 

To illustrate the model with non-zero reflectivities, the amplifier used in section () is simulated for 
two different set of reflectivities. In the first computation Rj and R 2 are both taken equal to 0. 1 . In 
the second one, Rj is taken as 0 and R 2 as 1 (double-pass amplifier). The electron density, photon 
density and point modal gain profiles corresponding to these amplifiers are illustrated in figure () 
and () respectively. 

From figure Q-i, It can be seen that the signal photon density is not continuous at z=0 (air-semicon- 
ductor interface). This is due to the denominator of equation (40) modeling the boost in the signal 
density inside the resonant cavity due to amplified internal reflections. Another discontinuity is 
present in signal photon density travelling in +z direction at z=L. This is just due to non-unity trans- 
mittance of the second facet. The signal wave in -z direction is weaker than the one in +z direction 
because for the given operational parameters, the one-pass gain of the device is smaller than R 2 
As a consequence, the electron density is more depleted towards the right hand side of the amplifier 
as seen in figure ()-ii. 

In the case of the double-pass amplifier, the photon density in +z direction reflect from the right 
facet of the diode (acting as a mirror) and amplified by the same one pass gain before exiting the 
amplifier from left end. Consequently the total photon density is higher toward the left facet of the 
diode depleting the electron population deeper in that region. 













Analysis of the effect of amplifier parameters on gain and efficiency 


The purpose of the modeling is to simulate the behavior of a quantum well laser amplifier as closely 
as possible in order to optimize amplifier parameters for maximum gain and efficiency and to un- 
derstand the trade-offs in trying to do so. 

The parameters mentioned above are length, width, facet reflectivities, confinement factor per well, 
number of quantum wells, excitation and pump power levels. 


Constraints on parameters involved in optimization 


a) Width 


To understand the effect of width on efficiency and gain of the amplifier it is necessary to look in 
equations ( 1 2), ( 1 3) and ( 14) more closely. Let’s assume that the series resistance of the laser diode 
is inversely proportional to its area, and RgQ is the series resistance in ohms of unit area.Then the 
series resistance R s is given by 



(18) 


Putting equation (12), (13), (18) in equation (14) yields 
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Now, if we assume that the active region of the amplifier receives the same input photon density- 
independent of the width of the amplifier, that is if the injected power is scaled proportionally to 
the width so that |a (+) | will remain constant, than the output photon density |c <+) | will remain 
unchanged. Consequently the efficiency will remain unaffected by the change in the width of the 
amplifier.This means that the problem of finding the optimum width and optimum input power can 
be considered together because the only parameters that is important is the power density injected 
to the amplifier. Thus, it is reasonable to keep the width constant and big enough to conform with 
the requirements of a broad area amplifier and try to find the optimum input power. A typical value 
for the width is 100|d .This value will be used through out the op timizati on, 

In reality, the electrical resistance of the diode is more accurate if a constant serial resistance is ad- 
ded to the right-hand side of (18) accounting for the effect of metallic contacts. This, results in a 
term in the denominator of equation (19) which is dependent on the width of the amplifier. The im- 
plication of such a change in the model is that the efficiency of the amplifier will decrease substan- 
tially with increasing width because of the increasing electrical power loss at metallic contacts. 

So, the width of the amplifier should conform with the requirements of a broad area amplifier but 
at the same time it should be kept as small as possible to increase the efficiency. 

b) Input Power 

High efficiency in a diode amplifier is achieved, if the pumped electrical power is effectively trans- 
formed into optical power. This requires operation of the device in a highly saturated regime, im- 
plying the use of relatively high excitation levels. However, the reflective facets of the amplifier 
can only operate up to a maximum value of power density beyond which they are damaged. This 
limit can be roughly set at (10 MW / cm^). Let’s denote this value by Dmax which stands for maxi- 
mum density. The corresponding maximum power limit P^^can then be expressed as 




In the case of a travelling wave amplifier the power density at the second facet is higher than the 
one at the first facet. Thus the input power level and the gain of the medium should be arranged in 
such a way that the output power density shouldn’t exceed the facet damage limit. 

The model doesn’t treat the spontaneous emission, the signal to noise ratio being very high for sa- 
turated operation regime below the threshold of oscillation. It has been reported [3] that if the in- 
put laser intensity is higher than 10 kW/cm , then regardless of the facet reflectivity, the output 
laser intensity is less dominated by ASE and less sensitive to the injected current density and input 
laser frequency. Accordingly, for a diode with a width of lOOjlm, confinement factor of 2% and 
quantum well thickness of 10 nm, the minimum injected optical power should be 0.005 W in order 
to neglect the photon density contributed by ASE. 


c) Confinement factor per well 


When an optical mode of the amplifier is too tightly confined to the active region, the divergence 
angle of the optical wave at the output facet is quite large resulting in considerable power loss in 
coupling light into the amplifier. Thus, there should be an upper limit for the confinement factor. 
Fallowing an empirical approach, this limit is experimentally determined to be 2.2%. 

However a small confinement factor is not very desirable either, due to decreasing efficiency resul- 
ting from a less saturated electron population. The excitation level can be increased up to some 
point to keep the input photon density to the amplifier constant, but the modal gain of the amplifier 
willstill be lower compared to the amplifier with higher confinement factor (1). 

The effect of other parameters will be explained in subsequent sections. 



Variation of gain and efficiency with pumped-in current 


The understanding of the changes in the performance of the device as a function of the excitation 
is very important because the understanding of the effect of the remaining variables can then be 
more easily achieved in this multivariate optimization problem. Also the excitation is continuous 
in nature and offers experimental ease in its manipulation compared to the injected power level, 
number of quantum wells etc. 

As it can be seen from previous graphs, the quantum-well gain rolls over at a high injected carrier 
density. This gives the amplifier the ability to operate at a high excitation level without exceeding 
the oscillation threshold. A positive consequence of this is low small-signal gain and high extrac- 
table power. 

The thickness of the active region of quantum-well laser amplifier being much thinner than that of 
a more conventional double heterostructure laser its internal efficiency is also considerably higher, 
resulting into a higher overall efficiency. A bigger portion of the electrical power being converted 
into optical power, the amount of heating is less severe, therefore contributing to a higher efficien- 
cy. The graphs below illustrates the variation of efficiency and net gain of the amplifier with in- 
creasing excitation. 

As can be seen from the first graph, amplifier efficiency initially increases with excitation. As the 
pumped-in electrical power increases, the fraction of input electrical power contributing to ampli- 
fier output optical power increases relative to the fraction spent in achieving transparency in the 
gain medium. Also, at low values of current density, the amplifier gain varies linearly, explaining 
the initial linear increase of the net gain with excitation in the second figure. The rate of increase 
of the efficiency decreases gradually with further increase in the current density. As a result, 






the efficiency goes up to a maximum /^(Peak efficiency) and then start going down. The main 
factor responsible from that change is the quadratic dependence of electrical power on current den- 
sity. With increasing excitation, the denominator of efficiency equation (electrical power) grows 
more rapidly than the numerator (extracted power). The gain roll-over and the saturation of the ma- 
terial gain makes the rate of increase of the extracted power with respect to current density (almost 
a constant for relatively high saturation) lower than that of input electrical power (quadratic). The 
third figure illustrates the variation of the electrical and extracted powers together with respect to 
current density. 

If the net gain at peak efficiency is not big enough, it can be increased by another 5dB at the ex- 
pense of decreasing efficiency (20%). This trade-off will be a part of our optimization algorithm 
letting us compare the gain of different amplifiers, their efficiency being decreased with increasing 
current density down to a minimum desirable efficiency (50%). 





Variation of gain and efficiency with injected power 


The effect of the saturation in the model can be better understood by looking at the variation of gain 
and efficiency for an amplifier with changing injection levels. The figures below illustrates the ef- 
fect of changing input power for an amplifier with relatively high pumping rate. 
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It can be observed that the efficiency increases very rapidly from 0 to 61 percent when the input 
power increases from 0.0001 to 0.6 Watts, where it reaches its peak and then starts going down with 
a relatively lower rate. The gain in the region corresponding to increasing efficiency is marked by 
an abrupt descent, the rate of decrease gradually going down with increasing injected power levels. 
It is interesting to note the asymmetry in the rate of variation of efficiency and net gain in the region 




of increasing efficiency. Eventually the explanation of that phenomena can be found in the numer- 
ator of the efficiency expression (the denominator being constant) by making the assumption that 
the diode acts like a travelling wave amplifier for reflectivities as low as 10^ (this assumption has 
been numerically verified). This is to say that the net gain of the amplifier is equal to one-pass gain 
Then we can rewrite the efficiency as, 


Eff = 


Pin(G s - 1) 



(23) 
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When the input power is very low, the gain of the system available to the injected photon is very 
high even in the presence of spontaneous emission due to non-saturated electron density. However, 
the small-signal gain is the maximum available gain and its value practically remain unchanged up 
to effectively saturating power levels. For a very small input power ( < lCf 6 W), despite the large 
value of the small-signal gain the product in the numerator of equation (23) is very small resulting 




into almost zero efficiency. This practically means that the extracted power is very low compared 
to provided electrical power although the gain is practically maximum. 

With increasing injected power, the electron population in the amplifier starts depleting according 
to equation (2), hence resulting into a lower net gain due to the logarithmic dependence of the gain 
on electron density. However, as far as the extracted power is concerned, the increase in the input 
power remains more important than the decrease of the net gain till peak efficiency is reached. 
Beyond that point, the decrease in the net gain is not overcome by the increase in injected power 
resulting into a gradual decrease in the efficiency. 

It can be observed that the decrease rate of the net gain with increasing injection level goes down 
The principal reason is the parallel decrease in the rate of increase of photon population resulting 
from increasing saturation. 



Variation of gain and e fficiency with the number of quantum wells 


The optimization of the number of quantum wells to be used in the amplifier is one of the main 
reasons why a numerical simulation of a detailed model was required. An experimental investiga- 
tion would be money-wise and time-wise very expensive. The numerical results summarized in the 
graphs below show that other input parameters being constant, there is an optimum number of 
quantum wells for which gain and efficiency are maximum assuming that the amplifier is mode- 
rately saturated. 
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As it is apparent from the first term in the right hand side of equation (2), increasing the number of 
quantum well decreases the unsaturated electron density per well. But in computing the total gain. 




we multiply the gain from each well by the number of quantum wells. And the effect of this linear 
multiplier is dominant over the effect of the reduction of electron population per well because the 
latter appears as a divider in the semilogarithmic gain expression. Thus we can conclude that the 
overall effect of increasing the number of wells is to increase the unsaturated gain. 

However, when the amplifier is operating at moderately saturated regime, as the graphs illustrate, 
an attempt to increase the gain will also increase the photon population present in the active region 
of the device, enhancing the saturation effect. Thus for sufficiently high number of quantum wells, 
the combined effect of the decrease in the first term and the increase in the second term of equation 
(2) reduces the electron population up to such extent that the linear multiplier (N^) in the gain 
expression is no more capable of increasing the net gain of the device. 
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The electrical and optical powers provided to the diodes are the same. The travelling wave ampli- 
fier approximation still being valid, equation (23) explains the similarity between the gain and ef- 
ficiency versus number of quantum well graph. 

For an undersaturated amplifier, increasing the number quantum well will result into an increase in 
both the efficiency and the gain contrary to the over-saturated case where an increase in the number 
of quantum well will result into a less efficient and lower gain device. 



Variation of gain and efficiency with amplifier length 


According to equation (10) the one pass gain of the amplifier (net gain in travelling wave amplifier 
case) is the integral of the of the net local gain along the amplifier. In the small signal case, the net 
gain of the amplifier increases linearly with respect to the length, the value of the local net gain 
being the proportionality constant. This simple treatment is no more valid for the large signal case 
due to the saturation effect. 


J-2500 A/cm-2 Pin=0.02W Eto=%2.2 NW=3 



The figure above illustrates the variation of the net gain of an amplifier with length for a moderately 
saturating input power and relatively high electron pumping rate. As it can be observed, the in- 
crease in the net gain of the amplifier with increasing length results in a higher photon population 
inside the active region enhancing the effect of saturation. This results in the roll over of the gain 




curve instead of the linear behavior predicted for the small signal case. However it is important to 
note that, beside the depletion of the electron population the net gain keeps increasing with length 
even for relatively long devices. 

To understand the behavior of the efficiency with changing length it is necessary to examine how 
the electrical and extracted powers are effected. The current density being kept constant, an incre- 
ment in the length results in an increment in the current provided to the diode. The series resistance 
of the amplifier however decreases due to the first expression in the right hand side of equation 
(14). 

After some simple manipulation the electrical power can be rewritten as 

P el = [(V f + R 0 J ) DJ] L+[R { (DJ ) '] L 2 (24) 

Where R 0 and Rj are the constants in equation (14). In the range of data the simulation is run, the 
first term in the right hand side of equation (24) is much bigger than the second one. This explains 
the linear relationship of electrical power with length in the graph illustrated below. 

The relationship established between the gain and the extracted power in equation (23), explains 
the variation of the latter with varying length. The second derivative of the extracted power curve 
is zero at a length close to 0.08 centimeters. At this point the second derivative is changing sign 
from positive to negative. The quotient of both powers (as indicated in the equation (15)) defining 
efficiency, reveals this phenomena producing a maximum at the same point which is nothing but 
the peak efficiency observed in the graph below. Thus, other parameter remaining constant, an am- 
plifier has an optimum length for which the efficiency is maximum 
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It is once again possible to increase the net gain at the expense of efficiency by increasing the length 
of the amplifier beyond its optimum efficiency value. 

The relative temperature variation corresponding to the analysis developed above is illustrated be 
low. The explanation of the shape of the curve is straightforward from equation (4). 
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Variation of gain and efficiency with optical confinement factor 


The optical confinement factor per well is one of the most important model parameters affecting 
the incident optical power density and modal gain at the same time. The parameter also appears in 
equation (2) but it cancels out with the one present in the modal gain expression (1). Eventually, 
it would be more convenient to use the material gain only to make more physical sense but it is not 
done so to preserve the notation used in Marcuse’s paper. 

To go from incident and exiting power levels to corresponding photon densities and vice versa, a 
classical approach is used relating the two by making use of the energy of a photon, the group ve- 
locity, the aperture of the active medium and confinement factor as expressed in equation (12) and 
(13). The use of confinement factor here accounts for the determination of the fraction of power in 
the active region of the quantum well structure, where amplified photonic energy will couple into 
the mode as described by equation (1). 

As far as the net gain of the amplifier is concerned the ratio of equation (13) by equation (12) will 
result in the cancellation out of the confinement factor present in both of the equations.The only 
parameter directly affected is then the one pass gain described by equation (10) which then mo- 
difies the signal photon densities at the diode boundaries as described by (8) and (9). 

It is straightforward from the corresponding expressions mentioned previously that increasing the 
confinement factor will result into an increase in the modal gain which in turn will increase the one 

pass gain. If the Fabry-Perot case, the signal photon densities at the boundaries of the device will 
\ 

then be higher resulting into a high overall photon population fallowing (5) and (6). However a 
higher photon density inside the amplifier will tend to deplete the electron population working 
against the increase of the net gain predicted from the previous analysis. The overall effect of in- 
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creasing the confinement can be seen from the graph displayed above which illustrates the roll over 
of the net gain following the saturation effect becoming more dominant over the counter influences 
mentioned before. 

It is interesting to point out that the net gain of the amplifier doesn’t increase substantially beyond 
2.2% confinement per well which was the limit set experimentally. However it is also important to 
note that both efficiency and gain can be simultaneously improved by forcing the previously set 
limit. Especially the gain in efficiency beyond the preset limit is far from being negligeable as it is 
clear from the graph presented below. 
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CONFINEMENT FACTOR PER WELL 


The behaviour of the efficiency can easily be understood recalling the travelling wave approxima- 
tion described by (23). 




Amplifier Optimization 


In the preceding sections, the variation of the gain and efficiency of the amplifier is described with 
respect to device parameters such as length, width, number of quantum wells, facet reflectivities, 
confinement factor and to operational parameters such as current density and injected power. The 
effects of the parameters are analyzed one at a time to give an intuitive feeling about the physics 
hidden behind the model equations. It is suggested that a travelling wave amplifier will perform 
better than a Fabry-Perot amplifier and that the confinement factor per well has to be as large as 
possible within the limit imposed by the coupling loss of the optical system. It is also shown that 
most of the parameters mentioned above has an optimum value maximizing efficiency others re- 
maining constant. Finally, some simple trade-offs in the optimization of efficiency and gain are out- 
lined. 

A simple optimization scheme would be to run the numerical simulation everywhere inside the 
physically significant region of the n dimensional parameter space, n being the number of param- 
eters to optimize. But even the most optimistic prediction of such task in terms of computation time 
is not less than a year. So an intelligent approach is needed to develop an empirical optimization 
scheme. One way of approaching the problem is to generate summary graphs which regroup the 
important information (peak efficiency) contained in a family of graphs. Such a method definitively 
causes loss of information but the cost of that is not very expensive if enough care is given in the 
analysis of the generated outputs. 

The procedure can be applied to two variables at a time. When the corresponding summary graph 
is constructed, the procedure can be extended to incorporate a third variable and so on 
The graphs below illustrate the variation of efficiency and gain with current density for amplifiers 
having different lengths. The first graph reveals an initial increase in the peak efficiency with in- 
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creasing device length fallowed by a decrease for relatively long amplifiers. The second graph por- 
trays the increase of one pass gain with increasing length. An increase in length results in a higher 
gain at all current densities. The gain curve seems to shift upward with a decreasing rate. This is as 
well due to the fact that with increasing length, equal length increments results in proportionnally 
smaller improvements in gain as to the saturation. The understanding of the effect of amplifier 
length on device performance would be simplified if extracted and electrical powers corresponding 
to the simulation presented above are analyzed. Previously mentioned graphs are presented below. 
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The fact that extracted power varies linearly with current density independent of the length (al- 
though with different slope) results from the fact that at corresponding injected power level the am- 
plifier is highly saturated and the increase in the photon population reduces an increase of output 
power down to a linear profile through depletion of electron population as mentioned in the previ- 




ous sections. In an amplifier operating in these conditions, an increase in the length will increase 
the gain of the device but with a decreasing rate as it can be observed from the slopes of extracted 
power curves. For very short amplifiers the net gain of the device is very small as it is clear from 
equation (10), no matter how big is the applied electron density in an experimentally reasonable 
range. 

The electrical power applied to the diode varies pseudo-linearly with current density for relatively 
short amplifiers as predicted by equation (24). It can be seen that the second term in the right hand 
side of the equation has a quadratic dependence on the length and thus is negligeable compared to 
the first term which only has a linear dependence. Furthermore, the forward voltage is much bigger 
than the product of the current density and R 0 even at very high current densities, resulting in an 
overall linear dependence, explaining the previous observation. With increasing 
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length, the terms quadratic in current density are no more negligeable leading to a parabolic vari- 
ation of applied electrical power with current density. 

For relatively short amplifiers, the rate of increase of injected electrical power catches up with that 
of extracted power at a higher current density, shifting the peak efficiency to the right. But with 
increasing length of the amplifier, the electrical power dramatically becomes a quadratic function 
of the current density, and thus forces the peak of the efficiency to smaller current densities. This, 
simultaneously explains the two-way variation of the peak efficiency in both current density and 
efficiency axes. 

It can be observed that the portions of the graphs corresponding to current densities below peak 
efficiency contains no useful information for the optimization because both gain and efficiency at- 
tain higher values at the remainding part of the simulation. Other portions of the graphs illustrates 
a previously mentioned trade-off between efficiency and gain. Even that information can be given 
away to only keep the value of peak efficiencies since a higher peak value will provide more room 
for a trade-off of efficiency with gain. The choice of quantum well laser amplifier instead of the 
more conventional double heterostructure laser was made fallowing a similar argument to increase 
the efficiency giving away a portion of the gain which would have been otherwise higher. Together 
with the peak efficiencies it is needed to keep track of the current densities and the gain correspond 
ing to these peaks to make sense out of the summarized information. The summary graphs gener 
ated are presented below. Remark that the peak efficiencies coincides with the small rectangles on 
the original efficiency graph. The simulation is run only up to peak efficiency thus saving a impor- 
tant amount of time. 

Now, it is possible to generate a family of summary curves varying a third parameter, injected po- 
wer or the number of quantum wells. As before, it is possible to shrink down the excessive informa- 
tion by only keeping track of the peak efficiency values of the graphs and the corresponding re- 
maining parameters. This procedure can be applied till the effect of all the variables is taken into 
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account. The summary of summary graphs obtained at the end should dictate many amplifier struc- 
tures with high gain and efficiency as desired. Note that the confinement factor has to be the last 
factor to be included since it may not necessarily result into a maximum in the summary graphs 
thus avoiding further shrinking. 

To give an example of the described procedure a sample overall simulation is implemented using 
the same gain function, width, transparency current density etc as before, but with varying number 
of quantum wells, confinement factor per well, amplifier length, input optical intensity and injected 
electrical power. 

To shorten the amount of simulation, the current density steps used to implement the first stage of 




summary graphs is kept relatively large resulting in plots with discontinuous derivatives at some 
points. This artifact is due to over or underestimating the current density at which the peak efficien- 
cy occurs and should be disregarded. 

The first page of figures show the variation of peak efficiency, peak gain, current density at peak 
efficiency and extracted and electrical powers at peak efficiency with amplifier length for different 
input optical power levels varying from 0.3W to 1.5W for an amplifier with a confinement factor 
of O.Ollper well and one quantum well in its active region. The second and third page illustrates 
the same simulation for devices with three and five number of quantum wells respectively. Fallow- 
ing that, the same analysis is repeated for a confinement factor per well of 0.016, this value being 
previously determined to be an experimental maximum. 

The explanation of these graphs is rather involved. Before this though, it is worth recalling that the 
optimization is aiming at predicting physical variables (length, confinement factor, width, and 
number of quantum wells) for amplifier structures that will yield optimum efficiency and gain at 
physically realizable input optical and electrical power levels. All of physical variables except am- 
plifier length requires very precise and irreversible epitaxial growth processes. Consequently, their 
prediction is more important than the others, for which the only requirement is to be within their 
physically realisable range. 

In figure plfl, it can be seen that increasing injected power increases the peak efficiency while 
shifting the length at peak efficiency to shorter values. The current density required at peak effi- 
ciency has become higher as well, while the peak gain has become lower. If one looks to (descrip- 
tion figure), it is clear that increasing the input optical power should require a higher current density 
to reach peak efficiency to produce the effect in figure plfl. This is in agreement with the satura- 
tion expression in the amplifier gain. 

This simulation is eventually more important than it seems. It shows that the most efficient ampli- 
fier is eventually very short for a given confinement factor. Increasing the injected power reduces 



the gain of the amplifier proportionnally. However, together with decreasing device length and in- 
creasing current density one can extract more optical power with respect to electrical power spent 
at the expense of gain. The material gain being very small for such a device and thus the overall 
gain being almost equal to unity, the photon and electron densities in the active region can be as- 
sumed constant along the amplifier. This result makes a lot of sense because it states that for an 
amplifier of infinitesimal length, there is an optimal set of electron and photon concentrations for 
which the amplification is maximum. However, such a device is far from being ideal since a gain 
of at least lOdb is required for an array made of these amplifiers to be of practical use due to cou- 
pling losses. So one has to increase the gain while keeping efficiency as high as possible. However 
this is not a very easy task because increasing the length of the amplifier will result into a nonuni- 
form electron and photon distribution along the device making the efficiency at an infinitesimal 
length element of the diode smaller than its optimum. The best set of parameters for optimum ef- 
ficiency for a given amplifier length is thus the one that on average makes the electron and photon 
concentrations along the amplifier deviate the least from their optimum value. Although, this pro- 
cedure is relatively easy to understand and visua;ize for a travelling wave amplifier, it is very dif- 
ficult for a device with non-zero facet reflectivities. The reason for that is the fact that in the latter 
case the photon and electron densities at a given location in the amplifier become dependent not 
only on the same quantities at previous locations but also on those in the remainder of the device 
length. That is the overall gain becomes a feedback factor as described by the Fabry-Perot model 
(see related equations). It is this complication that necessitated a simulation approach to this prob- 
lem in order to understand the physics of the device and to optimize its parameters, making an an- 
alytical treatment impossible. 

Returning back to optimization, one can now look at the remaining figures by keeping in mind that 
the purpose is to find amplifier parameters that results into a gain higher than lOdb while keeping 
the electron and photon densities close to their ideal value on the average. Figlpl and fig2pl shows 



the trade off between gain and efficiency very clearly. One has to remember however that more 
gain can be achieved by increasing either current density or amplifier length at the expense of ef- 
ficiency. This is where the choice of number of quantum wells and confinement factor comes in 
the picture. For a confinement factor per well of 0.011, an amplifier with 5 quantum wells and a 
length of 0. 135 centimetre seems to be an acceptable choice promessing an overall gain of 1 1 ,4db 
at 50% efficiency. This amplifier is predicted to achieve this performance when supplied with a 

9 

current density of 3200 A/cm and an input optical power of 0.3W. It is interesting to note that the 
lowest signal power offered the best trade-off between gain and efficiency yielding the highest ef- 
ficiency at the maximum gain with minimum electrical power. Also, as predicted previously higher 
confinement factor per well (0.016) yielded better gain and efficiency achieving 1 ldb at 53% effi- 
ciency with an input power of 0.3W and excitation of 2500 A/cm^ . It is clear that better results 
could be predicted by using lower optical power. 
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